間 分 割 多重 の 信号 処理 を FPGA に 実装 する 


スマ ー ト ・ ア ン テ ナ の 


ビー ム ・ フ ォ ー ミ ング 技術 


スマ ー ト ・ ア ン テ ナ 技術 は IEEE 802.11n で 注目 が 高まる 
MIMO (Multiple Input Multiple Output) 技術 と と も に , 
3GPP や 3GPP2, IEEE 802.16, IEEE 802.1 1 な ど , さ 
ま ざ ま な 通信 シス テム で 次 々 と 取り 入れ られ る と 考え られ る . 
本 章 で は , スマ ー ト ・ ア ン テ ナ ・ シ ステ ム の 明 本 原理 と 信号 処 
理 の 実装 に つい て 説明 する . (筆者 ) 


近年 , 音声 通話 以外 に 高速 な デー タ ・ サ ービス へ の 期待 
が 高まっ て いま す . プロ バイ ダ 側 で は , 基地 局 で 収容 可能 
な ユー ザ 数 を 増やし な が ら 設 備 費用 を 節減 し , 適切 な サー 
ビス を 開発 し な けれ ば な り ま せん . 
線 シス テム で デー タ 伝送 の 高速 化 を 実現 する に は , 制 
限 され た 周波 数 資源 の 利用 効率 を 向上 する 必要 が あり ます . 
この た め に 導入 され た の が スマ ー ト ・ アン テ ナ ・ シ ステ ム 
( Smart Antenna System) 和 1 で す . 


> 


注 1: アダ プティ ブ ・ ア レイ ・ アンテ Adgaptive Array Antenna) と も 呼ぶ . 


( a) 無 指向 性 アン テ ナ 凶 ( b) 120", 3 セク タ ・ ア ン テ ナ 凶 


1 従来 の 基地 局 アン テ ナ 


Minseok Kim 


1.。 ス マー ト ・ ア ン テ ナ と は 


人 @ スマ ー ト ・ ア ン テ ナ の 動作 イメ ー ジ 
割り 当て られ る 周波 数 帯 の 不足 に より , 制限 され た 周波 


数 帯域 内 で 通信 容量 を 増やす こと が 大 き な 課 題し と な っ て い 
ます 、 
図 て a) の よう に , 全 方 向 に 対し て 指向 性 を 持た な い ア 
ン テ ナ 基地 局 で は , 各 ユ ー ザ 信号 の 送受 信 は , 同一 周波 数 
を 利用 する セル 内 の ほか の ユー ザ に その まま 干渉 し て し ま 
いま す . 不要 な 方 向 に 無駄 な 電力 を 放射 する た め , その シ 
ステ ム 性 能 が 干渉 電力 に より 大 きく 制限 され る こと に な り 
ます . この 問題 を 解決 する た め に 工夫 し た の が セク タ ・ ア 
ン テ ナ で す . セク タ ・ ア ン テ ナ と は , 指向 性 を 持つ アン テ 
ナ を 用 いて カバ ー 領 域 を 分 割 し ます . 同一 セル 内 に お いて , 
異な る セク タ 間 の 干渉 を 抑え る こと が で き , 複数 の ユー ザ 
を 収容 で きる 基地 局 ア ン テ ナ ・ シ ステ ム で す . 図 【 b) は 3 
セク タ ・ ア ン テ ナ の 例 を 示し て いま す . 

スマ ー ト ・ ア ン テ ナ は , 空間 信号 処理 に より 同一 周波 数 
の 干渉 を 抑え る こと が で き , 制限 され た 周波 数 で シス テム 
容量 の 増大 を は か る 技術 で す . ある ユー ザ 方 向 に アン テ ナ 
指向 性 の メイ ン ・ ビ ー ム を 向け , ほか の ユー ザ に は ヌル を 
形成 し 送受 信 を 行う こと で , 干渉 の 影響 を 取り 除け ます . 
この 方 式 は 空間 的 に チャ ネル を 分 割 す る の で , 空間 分 割 多 
重 方 払 SDMA : Space Division Multiple Access) と 呼 
ば れ ま す . 図 2 に その 概要 を 示し ます . 

周波 数 の 利用 効率 を た 上 げ る た め に , 距離 の 離れ た 場所 異 
な る セル ) で は 周波 数 を 繰り 返し 利用 し て いま す . スマ ー 


ン テ ナ , 空間 分 割 多重 方 式 , 適応 ビー ム 形 成 , 
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図 2 

スマ ー ト ・ ア ン テ ナ 
を 用 いた 空間 分 割 多 
重 SDMA : Space 
Division Multiple 
Access) の 概念 図 


に 対し て 送 
お ける 干渉 電力 を 大 


ト ・ アンテナ 技術 を 用 いる と , ある 条 作 方 向 ) 
信 電 力 を 制御 で きる た め , 隣接 セル に 
幅 に 抑え られ ます . 

3 に , 携帯 電話 の セル ラ 基地 局 に お ける 伝搬 環境 の イ 
メー ジ を 示し ます . 伝搬 路 の 状況 に と より さま ざま な 現象 が 
起き ます . 建物 な ど に より 電波 が 反射 , 回 折 , 散乱 が 起き 
複数 の 伝搬 路 を 通過 する た め に 生じ る マル チ パ ス ・ フ ェ ー 
ジン Multipath Fadimg) や , 同一 周波 数 を 利用 し て いる 
隣接 の セル の ユー ザ に よる 同一 チャ ネル 和 干 涯 CCI: Co- 
伝搬 遅延 時 間 藻 遅延 広がり 特性) 
に 起因 する シン ボル 間 干 涯 ISI : 
に よっ て 通信 容量 や 品質 が 劣化 し て し まい ます . 

0 受信 アン テ ナ に 到来 する 信号 は , 図 人 4 a) の よ 

に マル チ パ ス に より , それ ぞ れ ある 位相 差 を 持っ て 辿り 
0 ます . その 信号 の 和 は 位相 差 に よっ て 強め 合っ た り , 
弱め 合っ た り し ます . この た め , 図 4 b) の よう に 受信 電 
力 が 激しく 変化 し , 時 間 あ る い は 場所 に より 通信 状況 が 劣 
化し ます . 最悪 は 途切れ て し まう こと に な り ま す . これ を 


channel Interference), 


Inter-symbol Interference) 


直接 波 較 


( a) マル チ パ ス 遅 延 波 の 影響 図 


スレ ッ シ ョ ルド 図 


フェ ー ジ ング 区 エー ジン グ 了 図 ェ エー ジン グ 罰 
区 間 図 区 間 較 区 間 較 


時 間 較 


図 4 マル チ パ ス 信 号 の 影 


基地 局 凶 


ョ こ 


遅延 広がり 


同一 チャ ネル 千 渉 図 と 


0 


マル チ パ ス ・ 


図 3 移動 通信 に お ける 伝搬 状況 


マル チ パ ス ・ フ ェ ー ジ ング と 呼び ます . 

マル チ パ ス ・ フ ェ ー ジ ング を 含め た 諸 問 題 に 対し て , 以 
前 か ら 送信 電力 制御 や 誤り 制御 符号 , ダ イ バ シ テ ィ 技術 が 
導入 され て いま し た . そし て , 周波 数 利用 効率 を より 高め 
る た め に , スマ ー ト ・ ア ン テ ナ を 用 いて 干渉 波 や 遅延 波 の 
除去 を する こと を 目的 と し た 研究 が 盛ん に 行わ れ ま し た. 
スマ ー ト ・ ア ン テ ナ の 適応 ビー ム 形 成 理 論 は , 古典 的 な 
適応 等 価 技術 と 等 価 で す . 最近 は , 送受 信 と も に 複数 の ア 
ン テ ナ を 考慮 し た MIMO シ ステ ム の 研究 に ウェ イト が 移り 
つつ あり , 一 層 活 発 に 研究 開発 が 行わ れ て いま す . 
スマ ー ト ・ アン テ ナ ・ シ ステ ム は , 図 5 の よう に アレ イ ・ 
アン テ ナ と 複数 の 無線 チャ ネル , 信号 処理 部 で 構成 され て 
いま す . ベー スバ ンド に お ける ディ ジタル 信号 処理 に より 
適応 的 な アン テ ナ の 指向 性 を 制御 で きま す . スマ ー ト ・ ア 
ン テ ナ に お ける ビー ム 形 成 は , 伝搬 路 の 状況 に 対し て 最適 
な 複素 重み 係数 を 計算 し , 各 ア ン テ ナ 素子 か ら の 信号 に 重 
み 係 数 を 乗算 し ます . これ は , 信号 の 振幅 と 位相 の 数 値 的 


ユー ザ 端 末 較 


Nmaa 


千 渉 信号 図 


マル チ パ ス 遅 延 信 号 較 


適応 等 価 図 


図 5 セル ラ ・ シ ステ ム の スマ ー ト ・ アンテナ の 動作 
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App 


App 


4 


な 調整 に な り ま す . この 処理 より , アレ イ 出 力 に は , 非常 
に 品質 の 高い K 最大 の MA: ノイ ズ と 干渉 電力 に 対す る 所 
望 信号 電力 の 比 ) 信号 が 得 ら れ ま す . 

適応 ビー ム 形 成 手法 と し て , 何 通 りか の 固定 ビー ム ・ パ 
ター ン を あら か じ め 用 意 し , 受信 信号 強度 に 基づい て ビー 
ム ・ パ ター ン を 切り 替え る 方 法 を , スイ ッ チ ド ・ ビ ー ム 
( Switched Beam) 形成 と いい ます . 一 方 , スマ ー ト ・ ア 
ン テ ナ に お ける 適応 ビー ム 形 成 Adaptive Beamforming) 
は , 伝搬 路 の 状況 に 応じ て 最適 な 重み を リア ル タ イ ム で 計 
算 し , 目的 の ユー ザ に 向け て メイ ン ・ ビ ー ム を , 干渉 と な 
る ユー ザ に 向け て ヌル を 形成 で きる の で , 信号 の 品質 が 改 
善 され ま す . 


@ 適応 ビー ム 形 成 に よる ヌル 走査 (ヌル ・ ス テア リン グ ) 

スマ ー ト ・ ア ン テ ナ は , アン テ ナ 指 向 性 パタ ー ン を 適応 
的 に 制御 で きる と 説明 し まし た . 特に ヌル の 走査 能力 が 大 
き な 特 徴 と な り ま す . ヌル を 走査 する こと で 干渉 流 の 影響 
を うま く 取 り 除け る か ら で す . 次 に , 指向 性 の パタ ー ン を 
制御 する 方 法 を 簡単 な 例 を 通し て 紹介 し まず すめ. 

簡単 の た め , 半 波 長 間隔 に 配置 され た 2 素子 アレ イ ・ ア 
ン テ ナ に 二 つ の 信号 が それ ぞ れ ある 角度 か ら 到来 する モデ 
ル を 考え まず 図 6). 仮に 信号 1 は 所 望 信 sj, 入射 角 : 
9 の) と し , 信号 2 は 干渉 信 景 sz, 入射 角 : 9。) と し ます . 
そう する と , ん 番目 の アン テ ナ に 受信 され る 信号 は 


こ 2z 

rr 570 em| テー サ ne | の) 
た 1 4 *《 1) 
2 


と な り ま す . ここ で , & り は 入射 波 の 包絡 線 , ん は 渡 長 , 
は 白色 ガウ ス ・ ノ イズ で す . さら に 簡単 の た め , ノイ ズ 


信号 1 


図 6 
2 素子 アレ イ と 2 波 の 信 
号 モ デル (⑰ 
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は な いこ と に し ます . アン テ ナ 間 隔 27 を ん 4/2 と する と , ア 
ゴイ ッ アン デ テ チ の 出力 は 。 


EEF ( 2) 
= 2 @ の 1 (7 2 eo2s(・exp(- 肥 sin6 ) 


と な り ま す . この 場合 , 所 望 信号 は 受信 し , 干渉 信号 に ヌ 
ル を 向け る た め に は , 以下 の 条件 を 満足 すれ ば よい こと が 
分 か り ま す . 


の 言 の zeKP( 7sin6 財 1 区 
の 時 の 2 交 *P(/rsin62F OK 


例え ば , 9 が 0' , 9s が - 30" の と き , 解 は , 


の =05-05/ 
の >=05+05/ 


と な り ま す . 

この と き , ビー ム ・ パ ター ン を プロ ッ ト し た の が 図 7 で 
す . ここ で - 830 方 向 に ヌル 落ち 込み ) が 形成 され て いる 
こと が 分 か り ま す . この ヌル の 方 向 を 制御 する こと が ス 
マー ト ・ ア ン テ ナ に お ける 適応 ビー ム 形 成 の 基本 原理 で す . 

メイ ン ・ ビ ー ム が 0" 方 向 か ら ずれ て いる の は , 最大 ビー 
ム 方 向 に 対し て は 制御 を 行っ て いな いた めで す . 最大 ビー 
ム 方 向 に つい て は , これ か ら 説明 する 平均 二乗 誤差 最小 化 
手 湊 MMSE: Minimizing Mean Square Error) に より 
最大 S7VA が 得 ら れる よう に 制御 で きま す . 


2. スマ ー ト ・ ア ン テ ナ に お ける 最適 人 


手法 


ここ で , スマ ー ト ・ アンテナ に お ける 適応 ビー ム 形 成 法 
と その 最適 化 手法 に つい て 説明 し ます . スマ ー ト ・ アン テ 
ナ の 機能 は , 適 店 ビーム 形成 と 適応 メル ・ ス テア リン グ 
( Adaptive NulLsteering ) が あり ます . 適応 ビー ム 形 成 は , 
受信 波 の 到来 方 向 が 未知 , ある い は 時 間 的 に 変化 する 場合 
に も , アレ イ の ビー ム を 自動 的 に 追従 させ る 機能 で す . 

一 方 , 強い 妨害 波 の 存在 下 で 微弱 な 所 望 電 波 を 受信 する 
場合 に 一 般 的 な 指向 性 合成 法 を 用 いる と すれ ば , 非常 に 低 
い サ イド ローブ ・ レ ベル を 設定 し な けれ ば な り ま せん . そ 
こ で , 指向 性 パタ ー ン の ヌル 点 を 自動 的 に 妨害 波 方 向 に 向 
ける 必要 性 が 生じ ます . これ が , 適応 メル ・ ス テア リン グ 


と いう 機能 で す . 
スマ ー ト ・ ア ン テ ナ は , 電流 環境 に 関す る 情報 を 学習 し 
な が ら , 指向 特性 お よび 周波 数 特性 を 環境 に 合わ せ て 変え 


て いき ます . 不要 波 と ノイ ズ を 含ん だ 電波 環境 か ら 所 望 の 
信号 を 抽出 する た め に , 所 望 の 信号 に 関す る 予備 知識 が 必 
要 で す . 通常 , 信号 の hb 財 波数 搬送 波 周波 数 ), 到来 波 
方 向 , 変調 方 式 , 偏 波 な ど が 利用 され ます . その た め , ス 
マー ト ・ ア ン テ ナ の 動作 原理 は , それ ら 予 備 知識 お よび 評 
価 関数 に よっ て 次 の よう に 大 別 さ れ ま ず 2. 
① 平 均 二乗 誤差 最小 化 洲 Minimum Mean Square Error: 
MMSE 
② 最 大 SNR 洲 Maximum Signalto-noise Ratio: MSN) 
③ 拘 束 付 出力 電力 最小 化 渡 Constrained Minimization of 
Power : CMP) 
④ 定 包絡 線 信 号 用 アル ゴリ ズム ( Constant Modulus 
Algorithm: CMA) 
ここ で は , 最も 広く 用 いら れ て いる 平均 二乗 誤差 最小 化 
注 MMSE 法 ) に つい て 簡単 に 説明 し ます . 


@ 平均 二乗 誤差 最小 化 法 (MMSE) 

MMSE 法 は , 所 望 の アレ イ 応 答 で ある 参照 信号 と 実際 の 
アレ イ 出 力 信号 と の 誤 六 誤差 信号 ) を 最小 に する こと に 
よっ て , 最適 な 重み 係数 を 決定 する 方 法 で す . この 方 法 は 
適応 メル ・ ス テア リン グ と 同時 に アダ プティ ブ ・ ビ ー ム 形 
成 を 行い ます . その た め に アン テ ナ 素 子 配 列 や アレ イ 構 成 
な ど に 制約 を 受け な いと いう 長所 が ある 一 方 , 参照 信号 と 
し て 厳密 に は 所 望 の 信号 その も の を 必要 と する と いう 矛盾 
が あり ます . 

実際 に は , 所 望 の 信号 の 性 質 周波 数 帯域 , 変調 方 式 な 
ど ) に 関す る 予備 知識 が ある の で , アレ イ 出 力 信号 を 適当 に 
処理 する こと に よっ て 適切 な 参照 信号 を 得る こと も 可能 で 
す . 従っ て , 受信 側 で 所 望 の 信号 の レプ リカ で ある 参照 信 
号 を 作る と いう 概念 は 現実 的 な 手段 で す . 完全 な 所 望 信号 
の 性 質 を あら か じ め 既 知 の 特性 と し て 解析 を 進め られ ます . 

最小 化 の 対象 と な る 誤差 信号 を 4 ), すなわち , 所 望 の ア 
レイ 応 徐 参照 信号 ) た の と 実際 の アレ イ 出 力 信号 * ) と の 
差 は 次 式 で 与え られ ます . 


e(= ァ r()-y( り = ァ 7( り りー @H (の SS Se ( 5) 


ここ で , gw は 信号 合成 に 用 いる 重み 係数 で す . 上 式 より , 
誤差 信号 の 2 乗 の 期待 値 平均 二乗 誤差 ) は 次 の よう に 表 さ 


れ ま す . 
呈 k の 用 = 中 の - 用 中 の -e"a の 用 @ 
MMSE 法 の 目的 は 重み ベク ト ル を 適切 に 選ぶ こと に 
よっ て 上 記 の 平均 二乗 誤差 を 最小 に する こと で す . 両辺 に 


重み ベク ト ル に 対す る 勾配 を 取り , 最小 点 を 求め る 問題 に 。 肖 
な り ま す . この 間 題 の 解 は 


の 


1 
二 全 折 た ABSN20028NBRSARSSSG3KeM59Ma (の 馴 聞 


と な り , 有名 な ウィ ー ナ 解 Wiener Solution) と 呼ば れ て 2 
いま す . ここ で , Au 7 は それ ぞ れ 信号 の 自己 相関 行列 , 
相関 ベク ト ルル 参照 信号 と 入力 信号 ベク トル と の 相関 ) に な 
り ま す . 


⑯ ディ ジタル 制御 アル ゴリ ズム 
MMSE 法 の 最適 な 重み 係数 が 以上 の よう に 与え られ る こ 募 

と が わか り ま し た . これ を 観測 可能 な デー タ か ら 求め る た 

め に , 近似 的 な 計算 を 用 いた 最適 化 ア ル ゴ リ ズム が 必要 と 

な り ま す . ディ ジタル 制御 に よる 最適 化 ア ル ゴ リ ズム と し 

て は , 最 急 降下 法 に 基づく LMS Least Mean Square) ア 

ル ゴ リ ズム と , サン プル 値 を 用 いた 直接 解法 の SMI 
( Sample Matrix Inversion) アルゴ リズム, さら に は 再帰 

的 最小 二乗 法 の RLS Recursive 中 0 アル ゴリ 

ズム が 広く 利用 され て いま す . 


ー80 -60 -40 -20 0 20 40 60 80 
角度 ”) 凶 


7 適応 ビー ム 形 成 に よる 指向 性 パタ ー ン の 一 例 
- 30" 方 向 に ヌル 落ち こみ ) が 形成 され て いる が, 所 望 方 大 0 方 向 ) 
に ビー ム が ずれ て いる こと が 確認 で きる 
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その 中 で , 還 確実 評価 関 


数 の 最小 点 に た どり つく こと , 計算 負荷 が 小さ いこ と な ど 
時 最 急降下 法 は 入射 濾 の 到来 角 が 接近 
し て いた り , 各 濾 の 電力 比 が 大 きか っ た りす る 場合 , 収束 
が 極端 に 遅く な る と いう 欠点 が ある こと が 知ら れ て いま す . 
この 問題 を 克服 する 方 法 の 一 つが SMI 方 式 で す . 

SMI 方 式 は 計算 負荷 が 大 きい こと か ら 従来 敬遠 され が ち 
で し た . 近年 , コン ピュ ー タ の 発達 に 伴っ て 注目 され る よ 
うに な り ま し た . RLS ア ル ゴ リ ズム は , 最 急降下 法 と SMI 
アル ゴリ ズム の 中 間 的 な アル ゴリ ズム で , 最適 な 値 に た ど 
りつ くま で の 時 間 収 東 時 間 ) が 短縮 で きる た め , 最も 一 般 
的 に 使わ れ て いま す . 


⑯ シス トリ ッ ク ・ ア レイ を 用 いる 実装 手法 

スマ ー ト ・ ア ン テ ナ の アル ゴリ ズム は , 回 路 実 装 に 適さ 
な い 逆 行列 や 割り 算 な どの 演算 が 必要 に な り ま す . また , 
高速 処理 を 実現 する に は , 演算 回 路 を 並列 化し た 効率 的 な 
構造 を 考え な けれ ば な り ま せん . この た め に 古く か ら 考 案 
され た の が シス トリ ッ ク ・ ア レイ Systolic Array ) 構造 で 
す . まず , シス トリ ッ ク ・ ア レイ に つい て 調べ て み ま す . 
1980 年 頃 , H.T. Kung と C.E. Leiserson は , VLSI 実 
現 向 き の ア ー キ テク チャ が 満た すべ き 条 件 と し て , 
① 少 な い 種 類 の 単純 な 構造 の セル を 用 いた 回 路 構成 
② 単 純 で 規則 正 し い デ ー タ と 制御 の 流れ 
③ 並 列 処理 , パイ プラ イン 処理 に よる 高速 計算 
③④ 入出 力 と 処理 の 重畳 
を 挙げ, これ ら を 満た す も の と し て , シス トリ ッ ク ・ ア レ 
イ を 提案 し まし た . シス トリ ッ パ ( Systolic) は , 心臓 収縮 
を 意味 する 形容 詞 で す . ちょ うど 血液 が 心臓 の 収縮 に 合 
せ て 体内 を 循環 する よう に , デー タ が クロ ッ ク に 同期 し て 
プロ セッ サ ・ ア レイ 上 を 移動 し て 処理 を 行い ます . 上 記 の 
条件 ① お よび ② は , 回 路 の 論理 設計 や レイ アウ ト 設計 , ア 


園生 園 中 国 


( a) 1 次 元 配列 較 


TTF ロロ 
(本 ロロ 


上 刻 貫 商 園 


図 8 

シス トリ ッ 

2 間 較 時 

の 構造 ( b) 2 次 元 正方 配列 図 ( c) 2 次 元 三角 配列 較 
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ル ゴ リ ズム 検証 を 容易 に する た め に 
張 性 の 点 か ら も 重要 で す . 

大 規模 な 回 路 を 一 つの チッ プ 上 に 実現 で き な い と き , 同 
ーー の チッ プ を 必要 個数 だ け 結 合 し て 一 つの 回 路 を 構成 で き 
る こと が 望ま し い 場 合 が あり ます . その 際 , 回 路 全 体 の 入 
出力 機構 が . その まま チッ プ 間 の 結合 に 使用 で きる こと が 
望ま し いで し ょ う . 

条件 ③ お よび 条件 ④ は , 高速 計算 を 実現 する た め に 重要 
で す . 特に , 処理 に 要する 時 間 が 比較 的 小さ な 場合 は , 
路 か ら の デー タ の 転送 と 処理 を 重畳 し , 入出 力 を 含め た 全 
体 の 処理 時 間 を 短く する こと が 重要 に な り ま す . 

シス トリ ッ ク ・ ア レイ は 多数 の 同一 ある い は 少な い 種 類 
の 簡単 な お プロ セッ シン グ ・ エ レ メ ン ト ( PE) を 規則 正しく 
接続 し て 構成 され ます . アレ イ の 形状 は , 図 8 に 示す よう 
に , 1 次 元 配列 , 2 2 次 元 三角 配列 が 基本 的 
な 構造 と な り ま す . タ は アレ イ へ 近 次 的 に 入出 力 さ れ , 
アレ イ 和 の 1 PE で 並列 に 処理 され ます . これ ら の 入 
出力 お よび 処理 は , パイ プラ イン 化 さ れ , 重畳 し て 進め ら 
れ ま す . 

シス トリ ッ ク ・ ア レイ お よび その 上 で の アル ゴリ ズム で 
ある シス トリ ッ ク ・ ア ル ゴ リ ズム は , 以下 の 条件 を 満た し 
ます 、。 

( 1) 均一 構造 

アレ イ は 同 -《 ある い は 少な い 種 類 ) の 簡単 な プロ セ シ ン 
グ ・ エレメント ( PE) か ら な り ま す . PE は セル ある い は シ 
スト リッ ク ・ セ ル と 呼ば れ , 各 セ ル は 有限 個 の レジ スタ を 
も ち ま す . 四則 演算 な か ど が 1 ステ ッ ズ 単位 時 間 ) で 実行 さ 
れ ま す . 

( 2) 局所 結合 と 局所 通信 

セル は 隣接 する セル と の み 直 接 接続 され まず 局所 結合 ). 
各 セ ル が 1 ステ ッ プ で 直接 通信 で きる の は , 近傍 の セル に 
限定 し , これ を 局所 通信 と 呼び ます . アレ イ 上 で の デー タ 
の 送受 信 は すべ て 局所 結合 を 通じ て 局所 的 に 行わ れ ま ず ほ 
か に , すべ て の セル に 共通 する 配線 と し て , 電源 と クロ ッ 
ク が ある が , これ ら は 局所 的 で は な い ). 

( 3) 逐次 的 入出 力 

初期 状態 に お いて , どの セル も あら か じ め デ ー タ を も っ て 
いま せん . 処理 され る デー タ は アレ イ へ あら か じ め 定 め ら れ 
た 一 定 の 割合 , 例え ば 1 ステ ッ プ に 1 個 あ る い は 2 ステ ッ プ 
に 1 個 ず つ 導 次 入力 され ます . この 割合 を 入力 パイ プラ イ 
ン ・ イ ンタ ー バ ル と 呼び ます . 出力 も 同様 に , 逐次 的 に 行わ 


重要 で す が , 回 路 の 拡 


回 


れ , 出力 の 時 間 間 隔 を 出力 パイ プラ イン ・ イ ンタ ー バ ル と 呼 


び ま す . この よう な 条件 を 満たす シス トリ ッ ク 構 造 は 
FPGA な どの ハー ド ウェ ア 実 装 に 適し て いる と いえ ます . 


3. FPGA を 用 いた 回 路 設 計 


ここ で は , FPGA を 用 いて MMSE 法 を シス トリ ッ ク 構 
造 で 実装 し た 例 を 紹介 し ます . 行列 の OX 分 触 QR 
Decomposition) を 用 いた 再帰 的 最小 二乗 法 の RLS ア ル ゴ 
リズ ム は , 高速 な 収束 特性 と 数 値 的 な 安定 性 が 知ら れ て お 
り , 今回 は この QRD-RLS を FPGA で 実装 する こと に し ま 
た 。 

9 に QRD-RLS シ スト リッ ク ・ ア レイ に 基づい た 受信 
シス テム の 概略 を 示し ます . RRK Radio Frequency) 帯 よ 

り 到来 し て きた 濾 は , アレ イ ・ ア ン テ ナ の 各 素 子 か ら 受 信 
され ます . DBRK Digital Beam Forming) 受信 機 に て IF 

( Imtermediate Frequency) 帯 の 信号 へ と 変換 され , AD 
変換 の 後 , 信号 は FPGA へ と 送ら れ ま す . FPGA で は 
DDG Digital Down Conversion) で さら に 複素 ベー スバ ン 
ド 信号 に 変換 され ます . QRD-RLS シ スト リッ ク ・ ア レイ ・ 
プロ セッ サ に よる 重み 係数 の 更新 が 行わ れ , 干渉 波 や 遅延 
波 の 影響 が 取り 除 か れ た アレ イ 出 力 を 得る と いう 仕組 み に 
な っ て いま ず ?. FPGA は 米国 Altera 社 の Stratb( EP1S40) 
を 使用 し まし た . 

複素 重み 係数 ベク ト ル は 各 ア ン テ ナ か ら の 受信 信号 に 掛 
けら れ ま す が , FPGA の DSP ブ ロッ ク を 用 いる と 高速 に 動 
作 す る 乗算 器 を 実装 で きま す . この よう な 特徴 を 生か し て, 


1 ゲイ ン 調 整 図 


A-D 凶 
A-D 図 


複素 乗算 器 を 実装 し ビー ム 形 成 処理 に か か る 配線 と ロジ ッ 
ク 回 路 を 低減 で きま す . 


@ 0 分解 型 RLS シス トリ ッ ク ・ ア レイ 

MMSE 法 の 制御 アル ゴリ ズム で ある RLS 法 は , 古く か 
ら , 等 化 処理 や 適応 フィ ル タ な ど , 数 多い 無線 応用 に 用 い 
られ て いま す . これ は , 一 般 に 次 式 の よう に 優 決定 
( Overdetermined : 式 が 未知 数 より 多い ) の 方 程 式 の 解 と 
し て 用 いら れ ま す . アレ イ ・ アンテナ の 場合 , アレ イ の 出 
力 は , 次 の よう な 関係 式 で 表現 され まず ヵ ン が). 


x(Ueo+xa(U の キー+ xy(0 の 。 ュ ミリ +e(1) 
x(2 の o+xe(2 の の +…+ xy(2 の コミ y(②+e(② て 9 


ゃ (が ) の o+ *2(7) の 1+… ry( が ) の ニモ y()e() 


ここ で , は 時 間 イ ン デ ックス , W は 素子 数 で す . こ 
は ベク トル 表記 に 書き 直す と , 


の よう に な り , & り の 平均 二乗 誤差 を 最小 に する 
求め る 問題 に な り ま ず !. 

図 10 に 44 素子 に お ける QRD-RLS シ スト リッ ク ・ ア 
レイ ・ プ ロ セ ッ サ を 示し ます . 丸い セル は 境界 セル , 正方 
形 の セル は 内 部 セル に 対応 し まず 図 11). 各行 で は , 境 
界 セ ル で 回 転 角 度 が 得 ら れ , その 角度 に 応じ て 内 部 セル で 
回 転 を 行う こと に よっ て の OX 分 解 が 行わ れ ま す . 最終 的 に 
各 セ ル に は , ヵ ) と の 値 が 更新 され 保存 され ます . 


li 
男 
Mo 
S 
RY 


図 9 

QRD-RLS シス トリ ッ ク ・ 
アレ イ ( 4 素子 アレ イ ・ ア ン 
テ ナ ) DDC: 陣 gital Down Conversion 


2M ア ッ _ 着 


QRD-RLS シ スト リッ ク ・ ア レイ 図 
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X1(5) xz(4) xs(3) Xa(2) y(1) 
x1(4) xz(3) xs(2) xa(1) 

%1(3) xz(2) xs(1) Y 

X1(2) xz(1) 


図 10 ー 
QR 分 解 を 用 い (⑥)ー- 6 


る シス トリ ッ 


・ ア ・ デ 
ロー ② ざ 一 


( 7) 図 


図 11 境界 セル と 内 部 セル 


重み 係数 ベク トル o は 次 に 示す よう な 後退 代 信 Back 
Substitution) と 呼ば れる 方 法 に よっ て 求め られ ます . 
QRD-RLS ア ル ゴ リ ズム の 流れ を 図 12 に 示し ます . 

各 セ ル は 図 13 に 示さ れる CPK CORDIC Processor 
Element) に よっ て 処理 が 行わ れ , また , その 処理 は 近 次 
形式 の CORDIC で 構成 され ます . 境界 セル と 内 部 セル は , 
CORDIC に お ける ベク トリ ング Vectoring Mode) と ロー 
テー ショ ン ・ モ ー ド (Rotation Mode) に それ ぞ れ 対応 し て 
いま す . 

ベク トリ ング ・ モ ー ド で は , 複素 数 x』 の 虚数 部 が の -CPE 
に よっ て 除去 され 実数 と な り , 次 に 9-CPE1 に よっ て xi 
が ゼロ と な り , 回 転 角 度 の の ぁ の , 9 が それ ぞ れ 求め られ ます . 
ロー テー ショ ン ・ モー ド で は , これ ら の 回 転 角 度 に 応じ て 
入力 ベク ト ルル RW x)」。 Tnm( )]『 を の -CPE に よっ て @ ほ 
ど 回 転 さ せ て , 次 に ヶ と x の 実数 部 と 虚数 部 を 9 -CPE1 と 
9 -CPE2 で 9 ほど それ ぞ れ 回 転 さ せま す 

パイ プラ イン 化 さ れ た 計算 プ ブロッグ CPE) は , デー タ を 
上 か ら 順に 流し 込む こと に より 効率 的 な 計算 が 行わ れ ま す . 

図 14 は , 4 素子 リニア ・ ア レイ ・ ア ン テ ナ に お いて , び 方 
向 に 所 望 波 が , - 40' 方 向 か ら 干渉 波 が 到来 する 場合 の 
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ーー 人 2 時 eat -* 


図 12 QR 分解 シス トリ ッ ク ・ ア レイ に お ける 重み 係数 の 計算 


Re x%) 図 Im 包 ) 図 


9 ECPE。 0 
し ec 


( a) ベク トリ ング ・ モ ー ド 図 
図 


Re x%) 図 Im 包 ) 図 


1 
6 較 In《 Xo7 ) 較 
( b) ロー テー ショ ン ・ モ ー ド 凶 
図 13 CPE( CORDIC processor element) の 構成 複素 形式 ) 


ビー ム ・ パ ター ン を 表し て いま す . 所 望 波 方 向 に メイ ン ・ 
ビー ム を 向け 干渉 波 方 向 に ヌル を 向け て いる こと が 確認 で 
きま す . 実装 し た RLS ア ル ゴ リ ズム が 干渉 波 の 影響 を 取り 
除く よう に 動作 し て いる こと が 分 か り ま す . 

リス ト 1 に , 2 素子 シス トリック ・ ア レイ の トッ プ ・ レ 
ベル VHDL 記述 を 示し ます . 図 10 を 参考 に 結線 を 確認 し 
て くだ さい . ベク トリ ング ・ ロ ー テ ーション ・ モ ー ド の コ 
ン ボ ポー ネン ト 記述 は 省略 し ます が , pp.106-111 の Appendix 
を 参考 に すれ ば , 容易 に 実装 で きる で し ょ う . 
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リス ト 1 2 素子 シス トリ ッ ク ・ ア レイ の VHDL 記述 例 
LTBRARY TEEE: 

USE TEEE.std 1ogio 1164 .a11: 

USE TEEE.sEtd 1og1o arith.ALL : 

USE TEEE.std 1og1o_UNSTGNED .a11 : 


ENTTTY sysEo11o 2ele TS 
GENERTC ( 
に 】 : Tnteder := 
) : 
PORT ( 
c1k : エ TN std 1ogio: 
rese 上 : エ N gtd 1ogio: 


X11 Re, X11 Tm 
X12 Re, X12 Tm 
Y1 Re, Y1 Tm 


エ N std 1og1iC_veotor(B-1 downto 
エ NY Std 1og1C Vecor(B-1 downto 
エ NY std 1og1C Vecor(B-1 downto 


R11, R22 OUT std 1og1o Vector(B-1 downto 0): 


R12 Re, R12 Tm : OUT sd 1og1C_ Vector(B-1 downto 
U1 Re, U1 rm 
U2_Re, U2 TITm 
) : 

END ENTTTY sysEo11o 2ele: 


OUT sd 1ogq1o veotor(B-1 downto 
OUT sd 1ogq1o veotor(B-1 downto 


ARCHTTECTURE エ ヒ ] OF sygto11oC 2e1e TS 


COMPONENT rotaEion mode 
GENERTC ( 
B : TnmEeder := 16 
) : 
PORT ( 
c1k : Std 1og1G: 
rese 上 : Std 1og1G: 


メ 1n H Std 1og1c veotor(B-1 downto 0): 
Y_in : エ N std 1og1o vector(B-1 downto 0): 
phi in エ N std 1ogio_veotor(B downto 0) : 

theta in エ N std 1ogio veotor(B downto 0): 


エ エ @ : OUT std 1ogio_veotor(B-1 downto 
ェ 1m : OUT std 1ogio_veotor(B-1 downto 
メ OuE 上 : OUT std 1ogio veotor(B-1 downto 
Y_Oou : OUT std 1ogiC veotor(B-1 downto 
ph ou 上 OUT std 1ogq1o Vector(B downto 
theta ou 上 OUT std 1ogic_ vector(B downto 0) 


) : 
END COMPONENT , 


COMPONENT Veotorind_mode 
GENERTC ( 


千 小 波 図 利 細 
- 40 
電 一 一 邊 ご 一 答 一 一 人 
図 14 ) り 20 ! 約 毅 
適応 ビー ム 形 成 倒 4 素子 アレ ーーー 。 
イィ . アン テ ナ ) | 272 アレ イ ・ ア ン テ ナ 図 


( 4) Yoshiaki Yokoyama, Minseok Kim and Hiroyuki Arai : 
7 の /e77e777277 の 77 の 7 Sys7o77c た S 44 の 7 が ve の 7 7 が 579 が PA g77 7 が 5 
Pe77o7yg72ce 契 vg7g7 が 7o7。。 2006 IEEE Vehicular Technology 
Conference( VTC 2006 fall), Montreal, CA, Sept. 2006. 


Minseok Kim 


東京 工業 大 学 


B : Tnteger 

) : 

PORT ( 
c1k : TN gtd 1og1o: 
rese ヒ : TN std 1ogio: 


x 1n : TN gd 1og1c vector(B-1 downto 0): 
Y_1n : T エ N gd 1og1c vector(B-1 downto 0): 


で ら : OUT std 1ogio_ veotor(B-1 downto 0): 
ph1 ou 上 OUT std 1og1o_ vector(B downto 0): 
theta ou OUT std 1og1c Vecor(B downto 0) 


) 
END COMPONENT : 


STGNAT, P12,O12 
STGNAT, P13 , P23,O13 , O23 
STGNAT, P14 , P24 ,O14 , O24 


td 1og1o_ Vector(B downto 
td 1og1o_ Vector(B downto 
gtQ 1og1ic_ Vector(B downto 


8TGNAL X22 Re,X22 Tm td 1og1o Veotor(B-1 downto 
8TGNAL Y2 Re,Y2 Tm,Y3 Re,Y3 1Tm 
gtd | 1ogio ・ 88EOB= 1 downto 


BEGTN 


BOUNDARY1 Yector1ng_mode 
GENERTC MAP(B => B) 
PORT MAP (cl1k, reset, X11 Re, X11 Tm, R11, PB12, 
012.) : 


TNTERNATL1 ェ oaion mode 
GENERTC MAP(B => B) 
PORT MAP (c1k, reset, XX12 Re, XX12 Tm, P12, O12, 
R12 Re, R12 Tm, X22 Re, X22 Tm, PB13, O13) : 


BOUNDARY2 YeCtorind_mode 
GENERTC MAP(B => B) 
PORT MAP (c1k, reset, X22 Re, X22 Tm, R22, P23, 
923) : 


TNTERNAL2 rotatton mode 
GENERTC MAP(B => B) 
PORT MAP (clk, reset, Y1 Re, Y1 Tm, P13, O13, 
U1 Re, U1 Tm, Y2 Re, Y2 Tm, BE14, Q14) 


TNTERNAL3 rotatton mode 
GENERTC MAP(B => B) 
PORT MAP (oc1k, reset, Y2 Re, Y2 Tm, P23, 023, 
U2_Re, U2 Tm, Y3 Re, Y3 Tm, BE24, Q24) 


囚 に 全 で 
っ ! 所 望 信号 較 


相対 感 廣 dB) 
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通信 回 路 に 役立つ 
CORDIC ア ルコ リ スム Minseok Kim 


CORDIC( Coordinate Rotation Digital Computer) ア ME 陸 | 間 
ル ゴ リ ズム は , 三角 関数 , 四則 演算 , 極座標 系 と 直交 座標 | 1 
系 の 変換 な どの 計算 の た め に , 米国 の Jack E. Volder 氏 に ! Cos@ー テ | 
より 1959 年 に 考案 され まし だ PND XY99 
CORDIC ア ル ゴ リ ズム は , ビッ ト ・ シ フト ( Bit Shift) と これ を も っ と 整理 し た の が 次 式 に な り ま す . 
足し 人 算 Add) の み で , ある 回転 角 度 に 対す る ベク ト ル 回 転 
を 反復 的 に 行え ます . 回 路 実装 が 容易 で ある こと か ら , さ ME 3245E2 《 ? 
ま ざ ま な 信号 処理 回 路 に 適用 され て いま す . 例え ば , ディ 
ジタル 検 濾 回 路 で の NCQ 数 値 制御 発信 器 ) や DDS Direct この 式 で cos の は tan ゅ で 変形 で きま す . CORDIC の 考 
Digital Synthesizer), 固有 値 / 特 異 値 分 解 , 各種 Lattice え 方 は , この tan ゅ を 離散 化す る こと に 着目 し ます . 要 す 
フィ ル タ , DFT( Discrete Fourier Transform) /DCT る に , 回 転 角度 6 を 離散 化し た 値 の 反復 回 転 を 行い, その 
( Discrete Cosine Transform) など, いろ いろ な 通信 回 路 回 転 方 向 に よる 符号 付き 和 で 表現 され ます . tan の は 以下 
に お いて 非常 に 役に立つ も の で す . ここ で は その 理論 を 解 の よう に 離散 化 さ れ ま す . 
説 し , VHDL で 回 路 を 設計 し て み ま す . 
に 《 3) 


の =tan~「 27 
1. CORDIC 法 の 原理 
これ は , 最初 の ベク ト ル か ら 45” 回 転 , 14′ 回 転 , 7 


図 1 に 示す よう に , CORDIC ア ル ゴ リ ズム は 基本 的 に あ 回 転 。… の よう 大 ま か な 回 転 か ら 段々 細か く ) 回転 方 向 
る 座標 Y ゴ xyj か ら " 革 デマ, ] に ある 角度 @ で 回 転 す を 制御 し な が ら ゅ に 近づい て いき ます . p ゥ ビット ・ ワ ー ド 
る 動作 を 反復 し ます . これ を 式 で 表現 する と , 回 転 の 公式 長 を も つ 入 力 ベ クト ルキ x,y] を , 角度 @ で 回 転 する 場 
を 使い 次 式 の よう に 書け ます . 合 , ヵ 1 回 の 反復 が 必要 に な り , 固定 小数 点 演算 の 場合 , 


ワー ド 長 で その 計算 精度 が 決ま り ま す . 
cos の を の =tan-《 2- め の 関 係 を 用 いる と 以下 の よう 


に な り , 

2 COS の ニニ 1 円 居 

4 JIRtan 絆 JR 821eieret8BNst| ( 4) 
"エッ ャ リッコ EE 5 に)o 
反復 式 に 書き 直す と 次 式 の よう に な り ま す . 
" キ xy] 2 科 
1 ee 

9 軸 の 遍 記 4 作 | ( 5 
図 1 ここ で , g ニ テキ 1( 回 転 方 向 を 示す ) と 巨 は 回 転 に より 生じ 


角度 の に対する ベク ト ル 回 転 る 利得 Processimg Gam) で あり , 次 の よう に な り ま す . 
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7 1 凶 
キキ キミ デニ 7 ブ cd 6 
11 二 放 0 
これ は , た に より 決定 され る 固定 の 値 な の で , ROM 


( Read Only Memory) など に あら か じ め 計 算 し て お く こ と 
が で きま す . また , それ ぞ れ の 回 転 に お いて ゲイ ン 調 整 を 
する こと で は な く , 最終 的 な ゲイ ン 過 で 調整 する こと で 計 
1 ます . 
2 に 反復 回 転 に よる ベク ト ルル 回転 の 様子 を 示し ます . 
2 回 転 方 ・ 向 は , 角度 アキ ュ ム レー タ 
人 語 の 関係 式 か ら 定義 され ます . 


キュ データー 呈 " tan~G の り 


( Accumulator) に よ 


ここ で , 角度 の 値 は LUI( Look-up Table) を 用 いて 実現 

され ます . 
で は , CORDIC ア ル ゴ リ ズム に お ける 二 つ の 基本 モー ド 

に つい て 調べ て み ま し ょ う . まず , Volder 氏 に より 名 づけ 
られ た ロー テー ショ ン ・ モ ー ド (Rotation Mode) と ベク ト 
リン Vectoring Mode) で す . ロー テー ショ ン ・ モー ド 
で は , 角度 アキ ュ ム レー タ を 回 転 角度 あと し て 初期 化し ま 
す . 各回 転 ス テッ プ に お いて , 角度 アキ ュ ム レー タ を 小さ 
くす る よう に 次 の ステ ッ プ で の 回 転 方 向 を 決め ます . こ 
よう な 動作 は 次 式 の よう に 表現 され ます . 


ロニ セー タ ・ 2 
5 に も / と ルン 
mm ニー ran_ (2 


( 茹 <0 の と きみ =-1, それ 以外 の と き g=+1) [ 


Se95K28S809875i ( 8) 


この よう な ステ ッ プ を ヵ 回 繰り 返す こと で , 以下 の よう 
な 角度 あの ベク ト ル 回 転 が 得 られ ます . 


ィ , 三 4( rocoszo 一 osn <o) 


坊 =4( yocOS<o 寺 roSim る o) 
x。 0 ( 9) 


4= Trz3 
7=0 


次 に , ベク トリ ング ・ モ ー ド の 場合 , 入力 ベク トル を ェ x 
軸 K y=0) に な る まで 回 転 し ます . この 結果 , x の 値 は べ 
クト ル の 大 き さ に な り , z の 値 は 入力 ベク ト ル の 位相 角 に 
な り ま す . ベク トリ ング ・ モ ー ド は , 各 ス テッ プ に お いて 
y の 値 を 小さ くす る よう に 次 の ステ ッ プ で の 回 転 方 向 を 決 
め ま す . ロー テー ショ ン ・ モ ー ド と 同様 に この 動作 は 次 式 


特 ワイ ヤレ ス 通 信 の 効率 を 高め る 


で 表現 され ます . 
エロ 1 ご ルー の 2" 
175 に ちよ 7 
L4 = ター 生 an (2 
( 叔 <0 の と き 図 =+1. それ 以外 の と き 図 = -1 ) 


SREEESSSHRS EE 《 10) 


また , 計算 ステ ッ プ を ヵ 回 繰り 返す こと で , 以下 の よう 
な 結果 が 得 ら れ ます . 


計 2 妥 
も 4 0) +)o 


=0 


サド 《 11) 
る 三 <o 填 tan「 (% の xo) 


4=] Tyrz? 
7=0 


ロー テー ショ ン ・ モ ー ド や ベク トリ ング ・ モ ー ド は 
初 の tan の 値 が 20 45” ) に な る こと か ら , - z/2 ぐ の @ ぐ ァ z 
/2 の 範囲 で 制限 され ます . 


2. 回 路 設 計 


CORDIC ア ル ゴ リ ズム は ビッ ト ・ シ フト と 足し 算 の み で 
実現 で き , 回 路 化 が 容易 で ある と 知ら れ て いま す . 
CORDIC ア ル ゴ リ ズム を 実装 する 方 法 は , 要求 性 能 や 回 路 
リソー ス 別 に いく つか 上 げ ら れ , シリ アル 加算 器 を 使う か 
と うか に より ビッ トド ・ シリアル / ビ ウト ~・ 
Serial/Bit-Parallel), また , 計算 ステ ッ プ を 回 路 で 反復 的 
に 再 利用 する か どう か に より 展開 / 反 復 Unrolled/Iterative) 

で 分 類 で きま す . 
に , 性 能 が も っ と も 優れ た Bit-Parallel Unroled 構 


図 2 ベク トル 回 転 を 繰り 返す 
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パラ レル Bit- 


図 3 Unrolled CORDIC 構造 


造 を 示し ます . 

CORDIC ア ル ゴ リ ズム で は , 利得 の 調整 スケ ー リ ング ) 
は 必ず 行わ な けれ ば な り ま せん . 前 述 の よう に ROM な ど 
を 用 いて あら か じ め 計 算 し て お ける た め , 計算 負荷 を 低減 


し 72 3 3 


4 k 番 目 の ダ ブル 回 転 ス テッ プ 


な り , さら に 固定 小数 点 に お ける 与え られ た 精度 で 次 の よ 
うに 近似 する こと で , 割り 算 の 実装 を 避け られ まず ?. 


74 


時 5 (の 
BIUE 4 放 Ci 《 19) 


この こと か ら , 利得 の スケ ー リ ング を ビッ ト ・ シ フト と 
足し 算 で 行え ます . 図 4 に 番目 の ダブ ル 回 転 ス テッ プ の 
ブロ ッ ク 図 を 示し ます . 

リス ト 1 と リス ト 2 で は , 回 路 の 効率 化 を 図る ダブ ル 回 
転 CORDIC ア ル ゴ リ ズム を 採用 し た VHDL 回 路 設計 例 を 


で きま す . し か し , その 値 を 用 いた 正規 化 に は 割り 算 や 平 
方 根 演算 を 必要 と する た め , 回 路 化 が 困難 と な り ま す . こ 
の 問題 は , ダブ ル 回 転 法 に より 解決 で きま す . 

ダブ ル 回 転 法 は , 角度 の @ で 回 転 す る の で は な く , 次 式 の 
よう に 角度 の /2 を 2 回 回 転 す る こと を 意味 し ます . 


ーー 2 
2 PR 《 12) 
gm ニタ 一 臣 an (2 
( 図 <0 の と き 了 図 =-1, やれ 以外 の と き 較 =+1 ) 
これ は , 各 ス テッ プ で 四 つ の ビッ ト ・ シ フト と 六 つ の 足 
し 算 で 実現 され ます . し か し , 利得 は 平方 根 な し の 状態 に 


Design Wave Mook 
動作 原理 , 設計 ・ 製 造 工程 か ら 応 用 事例 まで 


MEMS 開発 & 活 用 スタ ー ト アッ プ 


示し ます . 入出 力 は 16 ビ ッ ト と し て 設計 し まし た が , 利得 
を 考慮 し て 18 ビ ッ ト に し て あり ます . ロー テー ショ ン ・ 
モー ド と ベク トリ ング ・ モ ー ド ともに, 最後 の 部 分 で ス 
ケー リン グ を 行っ て いま す . 


参考 ・ 引用 \ 文献 
( 1) Volder J. E. : "77e COA//C か 7goo77e77C CO が 7 が 79 7eC7777 の 6 
IRE Trans action Electronic Computing, VolL EC-8, 1959. 
( 2) Jurgen Gotze, Steffen Paul and Matthias Sauer :" 4z 万 訟 
cze77 gco/7-/e 47go7777 767 Pg7g/7e/ 太 7gevg7e Co が の 727O7。 , 
IEEE Trans action on Computer, Vol. 42 No. 9, Sep. 1993 


Minseok Kim 
東京 工業 大 学 
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特 ワイヤ レス 通信 の 率 を 高め る 


リス ト 1 CORDIC ローテ ーション ・ モ ー ド の VHDL 記述 


ーー 天天 二天 天天 大 天天 天天 天天 天天 素 大 天天 天天 天天 誠 天天 天天 天天 天天 天天 ーー RETURN std 1ogio veoEor T8g 

CORDTC BTT-PARATLTET ROOTATTON VARTABTLE z:SEd 1od1C VecOr(x! エ ande) : 
BEGTN 
FEATURES TEF ds1gn='1! THEN -- 1 1g neqga ヒ 上 1ve 
12b1ts Aoouraoy と = ーー 
12bitgs エ /O ELSE -- 0O 1g pog1ive 
12bitsg atan LUT レン 本 す 】 
twice rota+on oF theta/2 for soa11ind w1thouE sqr ヒ 上 END TTF: 
RETURN (2) : 
theta 18 16Db1itg aoouraoy を END: 
p1 ーー テ END cord1c package 
p1/2 --> 0111 1111 1111 1111 ( 32767) 
0 --> 0000 0000 0000 0000 ( 0) LTBRARY TEEE: 
-p1/2 --> 1000 0000 0000 0000 (-32768) USE TEEE.std 1ogio 1164 .a11: 
USE TEEE.std 1oq1o ARTTH.a11 : 
d(rad) --> 1fF d>0, Fix(d/ (pi1/2)*( 2^(15-1) - 1 ) ) (16b+) USE TEEE.std 1ogq1ioc STGNED .a11 : 

1f d<0, fix(d/ (pi/2)*(-2^(15-1)  ) ) (16b) USE work.Cordic paokkage .ALL : 


, are 16-Di ヒ aoouraoy Du 18-Dbi word adqding 2bitg る ENTTTY oordio bp ro TS 

1 --> 00 0111 1111 1111 1111 ( 32767) GENERTC ( 

-1 --> 11 1000 0000 0000 0000 (-32768) B : Tnteqder := 16: 
PTPENO : natura] : 1) : 
x(k+1) (1-2^(-2k) ) *x(k) - d*2^(-k+1) *y(k) PORT ( 
マ (k+1) = (1-2^(-2k) ) *y(k) + d*2^(-k+1 ) *x (k) c1k : TN gtd 1og1o: 
Z (k+1 ) theta - dxaroctan(2^-k) ena : TN gd 1og1oj 


ーー Tnput zz : 2/2 : : TN gtd 1ogio Vector(B+1 downto 0): 
ーー 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 ーー 3 : TN gtd 1og1c vector(B+1 downto 0): 
LTBRARY 1eee: 送 : TN gd 1ogic vector(B-1 downto 0) : 
USE ieee.std 1ogio 1164.a11: 
x ou ヒ 上 : OUT std 1ogio_ veotor(B+1 downto 0): 
PACKAGE corQd1c package T8 Y_ou 上 : OUT gtd 1og1c vector (B+1 downto 0): 
FUNCTTON b1t]gsh1fF (xx: std 1og1C VeCor: 1:ma ヒ u エ a1 ) 2 oOu 上 : OUT std 1og1c Vector(B-1 downto 0 
RETURN sd 1od1C VeCor : ) : 
FUNCTTON b1trsh1F (xx: std 1Og1C VGCor: 1 :ma ヒ u エ a1 ) END ENTTTY Cord1io DD エ o: 
RETURN EQ 1o91C VectCr: 
FUNCTTON decis1gn (xx: std 1od1C_VeCtor: dgtgn: gtd 1og1o) 
RETURN sd 1od1C VGC : ARCHTTECTURE 上 ] OF oord1io bp ro TS 
END cord1io paolkage -- 14bitsg 
TYPE LUTARCTANTAB TS ARRAY(O to B-1) OF 
LTBRARY TEEE: gtQ 1og1C Vector(B-1 downto 0): 
USE TEEE.std 1ogio 1164 .a11 : CONSTANT ARCTANTAB : UTARCTANTAB := 
USE TEEE.std 1og1io ARTTH.a11: ( "01000000000000", "00100101110010", "00010011111101", 
USE TEEE.std 1og1o STGNED .a11 : "00001010001000", "00000101000101", "00000010100011", 
"00000001010001", "00000000101001", "00000000010100", 
PACKAGE BODY cordio paolkkage TS8 "00000000001010", "00000000000101", "00000000000011", 
FUNCTTON b1]sh1F (xx: std 1oOg1C VGCor: 1:ma ヒ u エ a1 ) "00000000000001", "00000000000001") : 
RETURN std 1og1o Vector TB8 
VARTABLE 2Z : Sd 1ogd1C Veo 上 or (x'! エ ande) : -- (1024, 604, 319, 162, 81, 41, 20, 10, 5, 3, 1, 1): 
BEGTN -- (16384, 9672, 5110, 2594, 1302, 652, 326, 163, 81, 
z(x'high) := x(x'high) : 作 1 。 20。 1 ユ 10 5。 3 
FOR 1 TN xx'high - 1 DOWNTO n LOOP 
Z(1) := (は - ロ ) : 8TGNAL xred, yred : Std 1og1c _ vector(B+1 dowmnto 0): 
END LOOB: STGNAT, zreg : gtd 1ogio vecor (B-1 dowmto 0): 
FOR 1 TN n-1 DOWNTO 0 LOOP 8TGNATL ssx, gy : Std 1ogio veotor (B+1 dowmnto 0): 
z(1) = !0!: 8TGNAT, 8z : Std 1ogio veotor(B-1 dowmnto 0): 
END LOOP: 
RETURN 2: BEGTN 
END : 


PROCESS8 (c1k, ena) T エ 8 
FUNCTTON b1trgsh1F (xx: sd 1Og1C VeCor: 1:ma ヒ u エ a1 ) BEGTN 
RETURN sd 1ogio veotor TS8 TF c1k'even モ AND clk = '1!' THEN 
VARTABLE z : Ed 1og1o_ Vector(x「 エ rande) : TEF ena = '1! THEN 
CONSTANT 1o: integer := x'high - n + 1: XX Out <= Xe: 
BEGTN Y_ OU 上 <= Y エ @edz 
TF 1o > 0 THEN Z OU <= 2Z エ @d: 
FOR 1 TN xx'!high DOWNTO 1o LOOP END TF: 
z(1) := xx(x'high) 』 END TF: 
END LOOP: END PROCES8 : 
FOR 1 TN 1o - 1 DOWNTO 0 LOOP 
(1) = (Tn) : メ (Kk+1 ) (1-2^(-2k) ) *x(k) - d*2*2^ 
END LOOP : マ (k+1 ) (1-2^(-2k) ) *y(k) + dQ*2*2^ 
ELSE Z(k+1 ) theta - dr*arotan (2^-k) 
FOR 1 TN xx'hiqgh DOWNTO 0 LOOP 
Z(1) := xx(x!high) : BSx <= deosiqn(b1trshif て (x in, PTPENO) , z inm(B-1 
END LOOP: By <= deos1gqn(b1trshifFt (y in, PTPENO) , z in(B-1 
END TF: Sg <= deos1ign (ARCTANTAB (PTPENO) , z in(B-1) ) : 
RETURN z: 
END : xred <= (( x in - DrghifE(x im, 2*PTPENO) ) - 
bit1ghifE (gy, 
FUNCTTON decis1gn (x: std 1og1C Vector: dgtgn: gtd 1og1C) Yreg <= (( y tn - bitrshift (y in, 2*PTPENO) ) 


) : 
) : 
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CORDIC ローテ ーション ・ モ ー ド の VHDL 記述 つづ き ) 


bit1sh1ft(sx, 1) ): 
Zred <= ZZ 1n - 82Z: 


END ARCHTTECTURE エ ヒ ]: 


TBRARY TEEE: 

USE TEEE.std 1ogio 1164 .a11: 
USE TEEE.std 1og1io ARTTH.a11 : 
USE TEEE.std 1og1o STGNED .a11 : 
USE work .Cordio paokage . ALL : 


ENTTTY cord1C _ rotation TS 
GENERTC ( 
B : Tntegde エ := 
NOofFPTPES : natu エ ra] := 


8 エ GNAL 


8 エ GNAL 


8 エ GNAL 


8 エ GNAL 


BEGTN 


x soale 1 3, x sca1e 7 9, x sca]1e 11 13 
Std 1og1o_ vector(B+1 downto 0): 
y_soale 1 3, y sca1e 7 9, y sca1e 11 13 
Std 1ogio veotor(B+1 downto 0): 


メ 8oa1e 
Y_soa1e 


CORLOOP : 
FOR k 1n 1 TO B GENERATE 
Cor : Cordio bp ro 上 


GENERTC MAP(B 


PORT MAP 


: Std 1ogic Vector(B+1 downto 0): 
: Std 1ogio_ Vector(B+1 downto 0): 


=> B, PTPENO => k-1) 


(c1k, ena, xx(k-1) , yy(k-1) , zz(k-1) , 


xx (k) , 


Yy(k) , zz(k) ) : 


: TN std 1ogio: 
: エ N std 1ogio: 
: エ N std 1ogio: 


: エ N sd 1og1o _ vector(B-1 downto 0) 


END GENERATE CORLOOP: 


xx(0) <= x 1n(B-1) g x in(B-1)& x in(B-1) & 
x in(B-2 downto 0): 
YYy(0) <= y in(B-1) g y in(B-1)g y in(B-1) & 
y_in(B-2 downto 0): 
zz(0) <= bitrshift(z in, 1): 


: エ N std 1og1o vector(B-1 downto 0): 
: エ N std 1og1o _ vector(B-1 downto 0): 


: OUT gtd 1og1c vector(B-1 downto 0): 
: OUT std 1og1icC_ Vector(B-1 downto 0) 


END ENTTTY Cord1o エ oa1on : 
ARCHTTECTURE エ て ヒ ] OF oord1o rotat1on TS 


COMPONENT cordio bp ro TS8 
GENERTC ( 
に 】 : Tnteder := 16』 
PTPENO : naEtura1 := 1)』 
PORT ( 
c1k : TN std 1og1oj 
: TN std 1ogio: 


: TN gd 1og1c vector(B+1 downto 0): 
: TN gd 1og1c Vector (B+1 downto 0): 
: エ N std 1og1io _ vector(B-1 downto 0): 


: OUT sd 1og1c vector(B+1 downto 0): 
: OUT sd 1og1c vector(B+1 downto 0): 
: OUT std 1ogiC_veotor(B-1 downto 0) 


END COMPONENT : 


TYPE TARRAY TS ARRAY(0 to B) OF stdQ 1og1C VeC モ or(B-1 
downto 0): 
8TGNATL, zz : TARRAY: 
TYPE OUTARRAY TS ARRAY(0 to B) OF std 1og1C_ VeCor (B+1 
downto 0): 
: OUTARRAY : 
: OUTARRAY : 


STGNAL, xx 
STGNAL yy 


-- gsoa11ng 
X_SCa 
X_8Ca 


X SCa 


X SCa 


Y_SCa 


x_ out 
Y_out 


e_11 13 


( 
( 
( 
( 


<= (x _ scale(B+1 ) 
<= (y_ scale (B+1 ) 


bitrsh1ifF て (xx(B) , 1) - 


bitrg 


b1trgh1ft (xx(B) , 7) - 


bitrg 


bitrgshift (xx(B) ,11) - 


bitrg 


btrgshif モ (yy(B) , 1) - 


bitrg 


bitrshifE (yy(B) , 7) - 


bitrg 


bitrshifE (yy(B) ,11) - 


bitrg 


x scale 1 3 - x soa]e 7 9) 
x soa1e 11 13 + bitrsh1ft (xx(B) ,15) ) : 
y_scale 1 3 - y soca]e 7 9) 
y_sca1e 11 13 + bitrshift (yy(B) ,15) ) : 


END ARCHTTECTURE エ ヒ ]: 


hift(xx(B) , 3): 
hift(xx(B) , 9): 


hift (xx(B) , 13) : 


hift (yy(B) , 3): 


hift (yy(B) , 9): 


hift (yy(B) ,13) : 


を x scale(B-2 downto 0) ) : 
を y Scale(B-2 downto 0) ) : 
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特 ワイヤ レス 通信 の 率 を 高め る 


リス ト 2 CORDIC ベ クト リン グ ・ モ ー ド の VHDL 記述 


ーー 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 大江 天天 大 ーー 


CORDTC BTT-PARALTLET,. VECTOR TNG 


LTBRARY TEEE: 

USE TEEE.sEtd 1ogio 1164 .a11: 
USE TEEE.sEtd 1og1o ARTTH.a11: 
USE TEEE.std 1og1c_STGNED .a11 
USE work .Cord1o paokaqde . AL : 


ENTTTY cordio bp veo TS 
GENERTC ( 
B : Tnteder := 16: 
PTPENO matu エ a] : 下 ) ば 
PORT ( 
c1k : TN sd 1ogio: 
ena : エ N gtd 1og1o: 


T エ N std 1og1o_Vector(B+1 downto 0): 
T エ N std 1og1o_Vector(B+1 downto 0): 
エ N std 1ogio Vector(B-1 downto 0): 


x_out 
y_out 
g ou 上 

) : 
END ENTTTY cord1io DD Veo: 


OUT std 1og1o Vector(B+1 downto 0): 
OUT std 1og1o Vector(B+1 downto 0): 
OUT std 1ogio vector(B-1 downto 0) 


ARCHTTECTURE エモ 上 ] OF Cordio bp veo TS 


-- 14big 
TYPE LUTARCTANTAB TS ARRAY(0 to B-1) OF 
Std 1og1C YeCor(B-1 downto 0): 
CONSTANT ARCTANTAB : UTARCTANTAB := 


2red <= 包 Tn - 82: 


END ARCHTTECTURE エ ヒ ]: 


LTBRARY TEEE: 

USE TEEE.std 1ogio 1164 .a11: 
USE TEEE.std 1og1o ARTTH.a11 : 
USE TEEE.std 1ogq1ioc STGNED .a11 : 
USE work .Cord1o paokkaqe . ALL : 


ENTTTY cordio arotan TS 

GENERTC ( 
B : Tntege エ 
NOoEPTPES maEura1 

) : 

PORT ( 
c1k : TN gd 1og1o: 
rese 上 : TN std 1ogio: 
ena : TN std 1ogio: 


x 1n : エ N gtd 1og1c vector(B-1 downto 0): 
Y_1n : TN gd 1og1c vector(B-1 downto 0): 


2_ou ヒ 
) : 
END ENTTTY Cord1C acan: 


OUT std 1og1o vector(B-1 downto 0) 


ARCHTTECTURE 上 ] OF oord1io arotan TS 


COMPONENT cordio bp Veo TS 
GENERTC ( 
に : TnEeqde エ := 16』 
PTPENO natura] := 1): 
PORT ( 


( "01000000000000『", 
"00001010001000", 
"00000001010001" 
"00000000001010", 
"00000000000001", 


(1024, 604, 319, 
(16384, 9672, 5110, 


8TGNAL xred, Y エ ed 


8TGNAL  zred 


STGNAL SXX, 8y 


STGNAL 


BEGTN 


"00100101110010", 
"00000101000101", 
"00000000101001", 
"00000000000101", 


"00000000000001"): 


162, 81, 41, 20, 


"00010011111101", 
"00000010100011", 
"00000000010100", 
"00000000000011", 


10, 5,。 3, 1, 1) 』 


2594, 1302, 652, 326, 163, 81, 


41, 20, 10, 5, 3, 1, 1) 


Std 1og1C Veoor (B+1 


Std 1og1c_ Vector (B-1 


Std 1ogd1C Veoor (B+1 


Std 1og1c_ Vector (B-1 


PROCES8 (c1k, ena) TS 


BEGTN 


TF oc]1k'event AND cl1k = '1' THEN 


エ EF ena = 


END TF: 
END TF: 
END PROCES8S: 


メ (k+1) (1-2^ 
Y(k+1 ) (1-2^ 


「1 ユ 1! THEN 
X Ou <= 
Y_Out <= 
2 Ou ヒ <= 


メ エ @d: 
YF@ed: 
g エ @G 


(-2k) ) *x(k) - dQ*2*2^ 
(-2k) ) *y(k) + d*2*2^ 


z(k+1) = theta - qrarotan (2^-k) 


Sx <= deois1gn(b1trsh1ft(x 1n, 
By <= deois1gn(b1trsh1ft(y in, 
8Z <= Qeoisgign (ARCTANTAB (PTPENO) , 
xreg <= (( zx in - bitrsh1ifE(x 1n, 

bi 上 1ghifE (By, 


yreg <= (( y in - bitrshift(y 1n, 
bi 上 1ghifF(gx 


(-k) *y (k) 
( -k) *x (k) 


PTPENO) , NOT 
Y_in(B-1) ) : 
PTPENO) , NOT 
Y_in(B-1) ) : 


NOT y in(B-1) ) : 
2*PTPENO) ) - 


2*PTPENO) ) + 


1) ) : 


1) ) : 


c1k : TN 
ena : TN 


TN 
TN 
TN 


x ou キ 
_Ou 
Z_ ou ヒ 
) 
END COMPONENT : 


std 1og1o: 
std 1og1o: 


Std 1og1C Vecor(B+1 downto 0): 
Std 1og1o_ vector(B+1 downto 0): 
std 1og1o veotor(B-1 downto 0): 


OUT std 1og1c vecor(B+1 downto 0): 
OUT std 1og1o vector (B+1 downto 0) : 
OUT std 1og1o vector(B-1 downto 0) 


TYPE TARRAY TS ARRAY(O to B) OF gtd 1og1C_ Veotor(B-1 


STGNAL zz 


downto 0): 
TARRAY : 


TYPE OUTARRAY TS ARRAY(0 to B) OF gtd 1og1c Vector (B+1 


STGNATL xx 
8TGNAL yy 


BEGTT 


CORLOOBP 
FOR k in 1 TO 


downto 0): 
OUTARRAY: 
OUTARRAY: 


B GENERATE 


COr : Cordio bp vec 


GENERTC MAP(B => 日 , 
PORT MAP (c]1k, ena, xx(k-1) , yy(k-1) , 


PTPENO => kk-1) 
ZZ (k-1) , 
メ (k) , yy(k) , zz(k) ) : 


END GENERATE CORLOOP : 


xx(0) <= xx 1n(B-1) g x in(B-1)g x in(B-1) & 


x 1n(B-2 downto 0): 


yy(0) <= y in(B-1) g y in(B-1) g y in(B-1) & 


Y_in(B-2 downto 0): 


zz(0) <= (OTHERS=>!0!): 


Z Out <= 2z2Z(B) 


ーー-bt1ghifFt (zz(B 1): 


END ARCHTTECTURE エ ヒ ]: 
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